Field induced superconducting phase in Superconductor-Normal 
metal and Superconductor-Superconductor bilayer 

X. MONTIEL and A. I. BuZDir^Q 
Condensed Matter Theory Group, LOMA, UMR 5798, 
Universite Bordeaux 1, 33405 Talence France 
(Date textdate) 

Abstract 

We study the proximity effect in a superconductor (S)-normal metal (N) bilayer systems under 
in-plane magnetic field and demonstrate that a compensation between the Zeeman effect and 
the energy splitting between bonding and anti-bonding levels may lead to a magnetic field induced 
superconducting phase well above the standard paramagnetic limit. It occurs that the non- uniform 
Fulde-Ferrell-Larkin-Ovchinikov superconducting state also exists in the field induced phase. The 
presence of the impurities scattering shrink the region of field induced superconductivity existence 
in S-N and S-S bilayers. 
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I. INTRODUCTION 



Quasi two dimensional superconductors have been studied for fifty years. A strong 
upper critical magnetic field (-^02) anisotropy was observed for the first time in interca- 
lated layered crystal of dichalcogenides of transition metal^. In these compounds, Hc2 
is higher for the in-plane orientation (ab plane) than in the perpendicular one (c-axis). 
Moreover, this Hc2 anisotropy was observed in intercalated graphite superconductors (in 
Cgi^ in CeCc^ and in Ce^fe^), in organic superconductor^, in high superconducting 
critical temperature (High T^) cuprates superconductord2'^i^ and also in superconducting- 
superconducting (S-S') YBa2Cu30T/DyBa2Cu30j and superconducting-insulating (S-I) 
(Y Ba2CuzOj)^ / {Vi Ba2CuzOj)^ artificial superlattice^^^"^. High Tc cuprates supercon- 
ductors have a layered crystal structur^^ and a strong electron anisotropj^^'^^'^^. The 
superconducting coherence length along the c-axis is smaller than the interlayer distance 
d. Consequently, high cuprate superconductors can be considered as natural superlat- 
tices. In high cuprate superconducting compounds, superconductivity exists in CUO2 
atomic planes which are sandwiched by non-superconducting atomic plane^^^'^. 

Ginzburg-Landau model (in the weak anisotropy limit (^c ^ (i))ff^and Lawrence-Doniach 
model (in the strong anisotropy limit (^c ^ give the description of the Hc2 anisotropic 
properties in layered superconductors near T^. This Hc2 anisotropy in superconducting 
multilayers can also be described microscopically by the standard Bardeen-Cooper-Schrieffer 
(BCS) and the tunneling Hamiltonian theory. Using this method, we obtain the {Hc2,T) 
phase diagram of layered superconducting systems. 

Some of high Tc cuprate superconductors can be considered as a stack of S-N, S-S^ 
or S-F^22l weakly coupled bilayers. The S-N, S-S' or S-F bilayer constitute the elemental 
unit cell of the multilayer. The properties of the S-N, S-S' or S-F bilayers qualitatively 
differs from a single S, N or F layers. Consequently, the properties of multilayers based on 
single layer elemental unit cell may be qualitatively different of the properties of multilayers 
based on bilayer elemental unit cell. We show in this paper that {H, T) phase diagram, 
with in plane magnetic field, of S-N and S-S' bilayers may reveal a magnetic field induced 
superconducting phase. 

The case of S-F multilayers has been studied iipSHlSl^ g_g i^jjayer, Buzdin et al have 
demonstratecP^"^ the possibility to overcome the paramagnetic limit at low temperature 
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FIG. 1: Mechanism of compensation of the Zeeman effect by the degeneracy hfting between the 
bonding and anti-bonding state. 
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FIG. 2: A superconducting (S)-Normal metal (N) bilayer with in-plane magnetic field H. 

for a high in-plane critical magnetic field. The field induced superconducting phase may 
appear at high magnetic field if the interlayer coupling energy t is higher than T^. Moreover, 
in this phase, the adjacent S layers have opposite signs of the order parameter (this is so- 
called 7r-stat^211)_ In this case, the Zeeman effect is compensated by the bonding- ant ibonding 
degeneracy lift produced by the hybridization between the two S layers (see figure [T|. The 
Cooper pairs in the tt state are more stable at high magnetic field than the 0-state. The 
0-state occurs when the adjacent S layers support the same signs of the order parameter. 
Somewhat similar idea in the context of two-band superconductivity was introduced by Kulic 
and Hofmann^-. In this paper, we show that in a S-N bilayer at high in-plane magnetic field 
H (see figure [2]), at low temperature and strong enough coupling t > Tco between the 
two planes, the paramagnetic limit is also enhanced above the usual Fulde-Ferrell-Larkin- 
Ovchinikov (FFLO)'^^''^ limit and a field induced superconducting phase may appear at high 
magnetic field. The corresponding mechanism is qualitatively the same as in the S-S bilayer 
but naturally there is no tt state realization in this case. We study also the influence of the 
impurity scattering on the {H, T) phase diagram of S-N and S-S bilayers. 

The outline of the paper is as following. In Sec. II, we present the model of a multilayer 
system and give the exact solutions of the Eilenberger equations. In Sec. Ill, we study the 



3 



influence of tlie transfer integral on the superconducting critical temperature and the effect 
of impurities. In section IV, we investigate the phase diagram of the S/N bilayer in the 
strong exchange field regime in both the clean and the dirty limits. In section V, we study 
the influence of impurities of the S-S bilayer (if, Tc) phase diagram. 

II. MODEL OF AN ATOMIC THICKNESS S/N BILAYER 

We start with a non-interacting model (see for exampld^^"^ of layered systems with 
alternating superconducting and normal metal layers. The electron motion is described in 
the N and S layers by the spin-dependent energy spectra ^n,a (k) and C,s,a (k) respectively. 
The parameters that characterize the systems are the transfer energy between the N and S 
layers t, the Cooper pairing constant A which is assumed to be nonzero in S layers only. The 
Zeeman energy splitting, due to in plane magnetic field H, is written as h = fisH where fis 
is the Bohr magneton. 

The two mechanisms destroying superconductivity under a magnetic field are the orbital 
and the paramagnetic effectP^ESl. Usually it is the orbital effect that is more restrictive. How- 
ever, in systems with a large effective mass of electron j^^**^^ * or in low-dimensional compounds, 
like quasi-one dimensional or layered superconductors under in-plane magnetic fielcP^, the 
orbital magnetism is weakened and it is the paramagnetic effect that is responsible for su- 
perconductivity destruction. 

The Chandrasekhar-Clogston paramagnetic limitP^"^ is achieved when the polarization 
energy of the normal electron gas, equals the superconducting condensation energy 

A^(0) Aq/2, where (0) is the density of state of the normal electron gas, Xn is its spin 
susceptibility, and Aq = 1.76Tc is the zero temperature superconducting gap. This criterion 
gives the exchange field hp{T = 0) = Ao/a/2 where the superconductor should undergo a 
first-order transition to the normal state. Larkin and Ovchinikov^ and Fulde and FerrelP^ 
(FFLO) predicted the existence of a non-uniform superconducting state with slightly higher 
critical field h^Ppj^QiT = = 0.755Ao > hp(T = 0). For quasi-2D superconductors the 
critical field of the FFLO state is even higher, namely h?pp^Q (T = 0) = Aq,'^ while in quasi 
one dimensional systems there is no paramagnetic limit at all^". We focus on the 2D case 
for which a generic temperature magnetic field phase diagram has been establishecP^. 

We consider the case when the coupling between the layers is realized via the transfer 
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energy t. In the whole paper, we assume t <^ Ep where Ep is the Fermi energy and 
then Cooper pairs are locahzed within each plane. The layers are coupled together by the 
coupling Hamiltonian 

Ht = tY. fc^i,^ (k) i^3,^ (k) + i^t (k) V^.+i,^ (k) + H.c] . (1) 

where ip^^ (k) (resp. ipj^a (k)) is the creation (resp. annihilation) operator of an electron 
with spin a and momentum k in the jth layer. In this paper, we study the S-N and S-S 
bilayers. In S-N system, the superconducting layer has the index j = and the normal 
metal j = 1. In the S-S bilayer, the superconducting layers are indexed j = and j = 1. 
The Hamiltonian of the system can be written as : 

H = Ho + HBcs + Ht, (2) 

where Ho is the kinetic and Zeeman Hamiltonian , Ht the tunneling Hamiltonian and Hbcs 
the BCS Hamiltionian. For the jth layer, the kinetic and Zeeman parts of the Hamiltonian 
are written as 

H^ = Y. fc-- (k, h,) (k) V^,,. (k)] , (3) 

(7,k 

The Zeeman effect manifests itself in breaking the spin degeneracy of the electronic energy 
levels according to 

0,. (k, h,) = (k) - ah, (4) 

where (k) = k^/2m — Ep i.e. for simplicity we choose the same electron spectrum in both 
layers. 

The field hj in the jth layer is assumed to be the same in both layers {h^ = hi = h) . 
We suppose an s-wave singlet superconductivity coupling which is treated in Hbcs within 
a mean field approximatiorP^' 

Hbcs = (k) ^,tt ("k) + A, (q) ^,,t (k) (-k)] + ^ / d^rA^^ (r) (5) 

where r is the two-dimensional coordinate within each layer and A the electron-electron 
coupling constant in the S layer only. The superconducting order parameter Aj is non zero 
only in the S layers as the coupling constant is in the N layer. In order to investigate the 
occurrence of modulated superconducting phase (FFLO), we choose the superconducting 
order parameter in the form 

A (r) = Ae^^-" 
5 



where q is the FFLO modulation wave vector. Using Gorkov's formahsm, we introduce the 
normal G and anomalous F Green functiond^ : 



(6) 



G,,, (k, k') = - (T. {ij^, (k) (kO) > = 5 (k - k' + q) G,, , (k) 
(k, k') = (T. . (k) ^+ (k')) > = 5 (k + k') (k) , 

where the brackets mean statistical averaging over grand-canonical distribution and the 
ordering operator in the Matsubara's formalism^^, and j and / the layer's indexes. From the 
equation of motionSH, the system of Green functions equation is in the Fourier representation 
in the S-N bilayer: 
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where u = {2n + 1) ttT are the fermionic Matsubara frequencies. In quasi 2D supercon- 
ductors, the maximal modulus of the FFLO wave vector is of the order of (^o)~^) being 
the typical superconducting coherence length. Since ^ ^ which of the order of the inter 
atomic distance with a good approximation we can consider (k + q) = ,^ (k) — h + \p.q^ 
where Vi? is the Fermi velocity vector in the plane. The anomalous Green function in the S 
layer writes 



-^0,0 



-A*A 



where A 



-aoA - /3f^ + t4 
{ioo -^{k) + h- vp.q) {ioo + ^ (k) + h), 
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{iu -^(k) + h- vp-q) {ioo + ^ (k) + h) and ^ = {iu - ^ (k) + h - vp.qf + iiuj + ^ (k) + hf . 

The superconducting order parameter in the 0th superconducting layer satisfies the self- 
consistency equation 



A* 

^0 



I^I^EE^oro = |A|T5^ / F+^e (7) 

To describe the FFLO modulated phase and the influence of the impurities it is more 
convenient to use the quasi-classical Eilenberger formalism. Moreover, we include the FFLO 
modulation phase and non-magnetic impurities. Applying Eilenberger's methocP^for layered 
systempS with Hamiltonian ([2]), the system of equations of Green functions can be written 
as: 
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where u = u + ih + (l/2r) and /^^^ (w, q) = ^ dC,F^i{u,^,q) is the anomalous 
Green function in the Eilenberger formahsm and r electron mean free pass time. We write 
ViT.q = vp.q. cos {(p) where (p is the polar angle (v^T^jq) and ()^ is the average over (j). We 
assume an in-plane scattering on impurities and the absence of spin flip during the electron- 
impurity interaction. To consider the presence of impurities we substitute u hj u + 1/2t 
and A* by A* + (fj^j (w, q))_^ /2r see for exampl^^l. 

Solving the Eilenberger equation ^ yields the Eilenberger Green function for the S layer 
labeled j = 

A* 



0,0 



2 1 



^)[f2if23+f22f23+2aia2] 
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(9) 



where we pose = + v'^q'^), fig = (a;^ + v'^q'^) with u = u + ih + (l/2r), uj± = 
u + ih + (l/2r) ±it, Us = CO — ivp-q- cos (0) with {ui^2 = W3 ± it = a3± — wp.q. cos (0)). The 
averaged solution on the (p angle of ^ writes 
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^) [01^3+^22^3+201^2] 
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(10) 



where ()^ is the average on the angle. Close to the superconducting critical temperature 
of the second order phase transition, the self consistency ([T]) can be writterP^ 



In 



T, 



CO 




(11) 



where Tc is the critical temperature of the superconducting layer in the S-N bilayer and 

2^2 

Tco = ^"^^^ e i^i^'^F the critical temperature of an isolated superconducting layer with m the 
electron's mass, kp the Fermi impulsion, 7 = 0.577215 is the Euler's constant and un the 
Debye frequency. 
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At zero temperature, close to the critical magnetic field of the second order phase ho, the 
order parameters are also small the self consistency ([T]) writes 

III. PROXIMITY EFFECT IN S-N BILAYER 

In this section, we investigate the superconducting phase in the S layer in the clean limit 
(r — oo). We study the superconducting critical temperature as a function of the interlayer 
coupling. We obtain the critical magnetic field of second order superconducting to normal 
metal phase transition as a function of the temperature and the interlayer coupling. Study 
of the influence of the impurities and of the S-S bilayers are respectively proposed in section 
IV and V. 



A. Critical temperature 

We study first the infiuence of the proximity effect on the superconducting critical tem- 
perature Tc of the S layer when no magnetic field is applied {h = q = 0) in the clean limit 



(r — )■ oo). Then (10) becomes : 



J 0,0 — 2a;(t2+a;2) ' K^"^) 



thus the self consistency equation writes 
In 



TA 1 



27 + 41n2 + ^ +^ - + 



2 27iTr \2 2tiT, 



where \l/ (x) is the digamma function. As seen in figjsj the superconducting critical temper- 
ature decreases with the increase of the proximity effect. At low transfer energy t <^ Tc, the 
superconducting critical temperature varies like = \ — |^ (t^) ' "'"^ ^^^^ 
interlayer coupling, the superconducting critical temperature reveals a quadratic decrease 
with the transfer energy. The superconducting state is not qualitatively influenced by the 
normal metal layer and can be considered as a single superconducting layer. 

At strong coupling between S and N layers aX t ^ Tc (but in the limit t <^ ujd) , the 
superconducting temperature varies as = The critical temperature decreases 

with the tunneling transfer as more and more Cooper pairs leak into the N layer. The 
superconducting properties in the N and S layers are practically the same and the bilayer 
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FIG. 3: Graph of Tc/Tco as a function of t/Tco (solid line). For t » Tco, the critical temperature 
of the superconductor decreases to zero. 

can be considered as an equivalent single S layer with an effective coupling constant A 
where A < A. In the case where t ^ wd, the S-N bilayer can be considered as a single 



B. Phase diagram of the S-N bilayer 

We study the {h, T) and {h, t) phase diagram of the S-N bilayer in the clean case and 

in presence of non-magnetic impurities. In a two-dimensional S monolayer, we can define 

three critical magnetic fields at zero temperature, = Aq/2 is the critical magnetic field 

for a second order phase transition, = A0/V2 is the critical magnetic field for a first 

order phase transition defined by Clogston-Chandrasekai'^^''^. h^^^^ = Aq is the critical 

magnetic field in the presence of FFLO modulations. One can see that h^^^^ > > ho. 

In a clean S monolayer with an applied in-plane magnetic field, the critical field is h^^^'~' 
[33] 



In this case, the Eilenberger anomalous Green function (|9]) becomes for arbitrary inter- 
layer coupling t : 



superconducting layer S with A = ^ as predicted irpH. 
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(14) 



where we note the appearance of three energy scales E^^ = h+lf'p.l^ and E12 = h±t+lt p.l^ . 
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1. {h,t) phase diagram at zero temperature 



From (14)and the self consistency equation (12) , the critical magnetic field h is shown 
to satisfy 

2 



hr 



t+ \l\{hc-tf -{q-vpf 



hc + t + 



\J\{hc + tf - {q.VFf\ . he + \J\hl- iq-vp) 



(15) 



where one must find the value of q that maximizes the critical field he- If the field induced 



phase is assumed to be uniform in the each planes, namely if g = 0, equation (15) merely 
reduces to 

|2 ^1 , ^1 ,4 



\hc\ . \hc — t\ . \hc + t\ = h 







(16) 



The number of solutions with physical meaning of the equations (16) differs with the value 



of t (see figure [4]) . We defines the critical interlayer coupling = V^ho = 1.2473Tco that 
determines the number of physical solutions. 



If t < tc, the equation (|l6|) has only one solution. The critical magnetic field at zero 
temperature writes hd = 



'2t2 + 2A/t^T4/;^. In the limit t < T^o, the solution can be 
written hd = ^ ^1 + ^ j • We note that the critical magnetic field at T = OK in the S-N 
bilayer increases with the interlayer coupling t. 



In the case t > t^, the equation (16) has three solutions with physical meaning. The 



first solution is hd- The second and the third solution are h 



c2 



2t2 + 2^yt^ - Ahl and 



/ic3 = \\jlt'^ — 2^/W-^4h^ respectively. In the limit t ^ T^o, the three solutions can be 
written as hd,2 = t ± AQ/32t^ and h^s = Aq/M. In the limit t ^ T^q, is of order of T^o/t 
and then h^s is of order of T^. Consequently, hc3 define the lowest critical magnetic field. 
For t = tc the critical fields hc2 and hcs coincide. 

In the case of high interlayer coupling t > tc , a. field induced superconducting phase 
appears at high magnetic field. This phase exists between the two magnetic fields hd,2 = 
t± AQ/32t^. Thus, the new zero temperature paramagnetic limit hd = t + A (0)^ /32t^ may 
be tuned far above the usual one h^^^^ = Aq merely by increasing the interlayer coupling. 



Thorough analysis of equation (15) shows that the upper critical field is even increased 



by an in-plane modulation (see figure [s]). 

The FFLO paramagnetic limit of the S-N bilayer also depends on the interlayer coupling 
t as seen in the figure [5] The field induced superconducting phase is observable at T = OK, 
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FIG. 4: (h/Tco, t/Tco) diagram for the S-N bilayer in the clean Umit (r — )• oo)at T=OK (sohd hne). 
The uniform superconducting state is presented in grey. The hne h = Ao/2 presents the critical 
magnetic field for a second order superconducting phase transition for a single superconducting 
layer. The line h = Aq corresponds the FFLO paramagnetic limit for a single superconducting 
layer. The line h = Ao/\/2 represents the first order paramagnetic limit for a single superconduct- 
ing layer. 

only when hc2 and hcs are distinguishable. In presence of FFLO modulation, the critical 
magnetic field at zero temperature h^-f^^ and h^^^^ are separated in the case t > I.STcq. 
Below this value, the usual superconducting {h, T) phase diagram may be strongly deformed 
(see figure [?]) . 



2. {h, T) phase diagram 



In this section, we study the second order {h, T) phase transition diagram taking into 



account FFLO modulation. The self consistency equation (|ll|) is 



Inff^ 

p27r 



Re 



Re 



7 + 21n(2) + |vl> i + 



i{h+vp .q. cos((/))) 
2^ 



+ 



i{h+t+vp.q. cos{<p)) 



I* u 



i{h—t+vp.q. cos((/))) 



(17) 



This analysis in general case can be performed only numerically on the basis of the equation 



(17). 



A magnetic field induced superconducting state appears at high magnetic field as we can 
see in the figure [6] for t = 2.Tco and |8] for t = S.Tcq. For /i ~ t, the Zeeman effect that 
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FIG. 5: {h/Tco,t/Tco) diagram for the S-N bilayer in the clean Hmit (r — )• oo) (solid line). The 
uniform superconducting state is presented in the grey region. The non-uniform superconducting 
(FFLO) phase in the S-N bilayer is presented in the dotted region. The line h = Aq/2 presents 
the critical magnetic field for a second order superconducting phase transition for a single super- 
conducting layer. The line h = Aq represents the FFLO paramagnetic limit for a single super- 
conducting layer. The line h = represents the first order paramagnetic limit for a single 
superconducting layer. 



tyTco=2, Second order line 
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FIG. 6: {h/Tco,Tc/Tco) phase transition diagram calculated for t = 2Tcq with the second order 
transition line (solid line) and FFLO state to normal state transition line (doted line). The inset 
presents a zoom of the superconducting re-entrance phase around /i ~ t ~ 'i.TcQ. 
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FIG. 7: {h/Tco,Tc/Tco) phase transition diagram calculated for t = 1.35Tco with the second order 
transition line (solid line) and FFLO state to normal state transition line (doted line). We see 
below Tc — 0.2Tco that the transition line is deformed. The compensation between the Zeeman 
effect and the bonding and antibonding state becomes relevant at low temperature. 




0,00 0,05 0,10 0,15 0,20 0,25 0,30 

Tc/Tco 

FIG. 8: {h/Tco,Tc/Tco) phase transition diagram calculated for t = 3Tco with the second order 
transition line (solid line) and FFLO state to normal state transition line (doted line). The inset 
presents a zoom of the superconducting re-entrance phase around /i ~ t ~ 3Tco. 

destroys the superconductivity is compensated by the bonding-antibonding states degen- 
eracy created by the proximity effect between the S and the N layers (see figure [T]). The 
lower and upper critical lines merge at field h = t and the field induced superconductivity 
is confined to temperature lower than Tm = vre^'^Tj)/ (8t) in the limit t ^ Tcq. Therefore, 
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the superconducting field induced phase is confined to temperature lower than Tm- These 
results were obtained for relatively strong coupling. For lower coupling, {t ~ Tco), the usual 
phase transition diagram is strongly deformed as shown in figure [7] and finally disappear for 
t smaller enough than Tcq. From an experimental point of view, one might choose a system 
with an intermediate coupling t small enough to settle superconducting field induced phase 
but large enough to separate re-entrance and usual S phase. 



IV. EFFECT OF THE IMPURITIES ON THE FIELD INDUCED SUPERCON- 
DUCTING PHASE 

In this section, we investigate phases with uniform superconductivity in the S layer. We 
study the superconducting critical temperature as a function of the interlayer coupling. We 
obtain the critical magnetic field of second order superconducting to normal metal phase 
transition as a function of the temperature and the interlayer coupling. Study of the infiuence 
of the impurities and of the S-S bilayers are respectively proposed in section IV and V. 



A. Critical temperature 



We start with the analysis of the infiuence of the impurities on the superconducting 



critical temperature. Then (10) writes 



/cto 



A* i +2 



(18) 



in accordance with the model developed irpH.The self consistency equation (11) in this case 
is written as 



In 



2.T y i^^^^iMl. 



(19) 



In the case of weak proximity effect {t <^ Tco), the decrease of the critical temperature is 



deduced from the equation (19)and reads 



In the clean limit {Tcq ^ t ^ ^) the superconducting critical temperature varies as 



1 



7 C(3) _ jT 1_\ 

8 7r2 192 tT,o J T'^o ' 



and the impurity scattering inside the N layer decreases the prox- 



imity effect. In the dirty regime (T^o ^ ^ ^ t) the superconducting critical temperature 
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FIG. 9: Graph of T^/Tco as a function of t/Tco- The clean case (l/2r = 0) is presented by the 
sohd hne. The impurity are plotted with repectively (1/2t) /Tcq = 0.01( 0.05, 0.097, 1 in dotted 
(respectively dashed, dashed-dotted and dashed-dottted-dotted) line. We see that the impurities 
enhance the superconducting transition temperature for weak interlayer coupling. On the other 
hand superconducting critical temperature decreases quickly in the presence of impurities at strong 
interlayer coupling. 

varies as = ^ — ^j^- The presence of impurities enhanced the superconducting state 



and Tc decreases slower than in the clean case (see figure [9]). In this case, the impurities 
decreases the effective transfer coupling and then the proximity effect. 

However at strong interlayer coupling t ^ Tco and l/2r ^ T^, the expression for the 
anomalous Green function ^ becomes /q+q = {} + (^M^^^^^wA^^ (^)) 

the critical temperature varies as = ^'^-^ ^1 — l^l^j- This means that scattering on 
impurities strongly decreases for high interlayer coupling as seen in the figure |9} In the 
regime t ^ Tco, the mixing between the superconducting state in the S layer and the normal 
state in the N layer is very strong. The bilayer draws near the regime A — > A/2 where the 
S-N bilayer can be considered as a single S layer with an effective coupling constant A < A. 
Note that Tc depends on the impurities contrary to the Anderson theorem prediction which 
is not astonishing because the system is non uniform. 
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B. Effect of the impurities on the phase diagram. 



In this section, we study the influence of the impurities on the [h, T) and (/i, t) phase 
diagram of the S-N bilayer. In the presence of impurities, the modulated phase disappears 
and h^^^^ decreases to h^^^^. When the normal phase is overcooled then the critical mag- 
netic field decreases from to h^. For simplicity in the whole paper, we will focus on the 
second order transition critical field of the S-N bilayer, taking in mind that if the transition 
is of order the first order the calculated field corresponds to the overcooling field and the 
critical region of superconductivity phase existence may be somewhat larger. Consequently, 
we study the influence of the impurities in the homogeneous case (g = 0). In this case, the 



anomalous Green function is the same as (18) with the substitution uj — oj + ih and can 
be written as 



1. {h,t) phase diagram at zero temperature 

The impurities change the form of the {h, t) phase diagram at T = OK as shown pre- 



sented in the figure 10 The {h, t) phase diagram has been calculated numerically. The 
critical interlayer coupling tc increases with the impurities diffusion potential l/2r. The 
maximal values of hd and hc2 decreases with the impurities diffusion potential contrary to 
/ic3- The variations of /id, /ic2 and /ics reveals that the superconducting phase in the S layer 
is enhanced by the presence of the impurities whereas the field induced superconducting 
phase is destroyed by the impurities. 

2. {h, T) phase diagram 

The (/i, T) phase diagram as been calculated numerically. The reentrance phase is 



strongly influenced by the presence of the impurities as seen on the figure 11 The maximal 
critical temperature under which the field induced phase exists, decreases with the impurity 
scattering potential. Moreover the upper and lower critical fields of the re-entrant supercon- 
ducting phase [hci and hc2) become closer with the increase of impurities diffusion potential 
as seen in the last part. In the case t = 2Tco, the reentrance phase totally disappears for an 



impurity diffusion potential l/2r upper than 0.097Tco. In the figure 11, the usual supercon- 
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FIG. 10: {h/Tco, t/Tco) diagram for the S-N bilayer in the clean case (r — )■ oo) (solid line) and with 
respectively (l/2r) /Tco = 0.05 (0.097, 1) in dashed (respectively dotted, dashed-dotted) line. 



(1/2t)/Tco=0 




(1/2t)/Tco=0,01 
(1/2t)/Tco=0,05 
(1/2t)/T«=0,097 
(i;2T)/Tai=1 



0,0 0,1 0,2 0,3 0.4 0,5 0,6 0,7 

Tc/TcO=0 



FIG. 11: (/i/Tco, Tc/Tco) phase transition diagram for the S-N bilayer calculated for t = 2Tco with 
the second order transition line in the clean case ((l/2r) /Tco = 0) (solid line) and with respectively 
(l/2r) /Tco = 0.01(0.05, 0.097, 1) in dotted (respectively dashed, dashed-dotted, dashed-dotted- 
dotted) line. The inset presents a zoom of the superconducting re-entrance phase around /i ~ t ~ 

ducting phase is also influenced by the presence of impurities. The critical magnetic field at 
zero temperature hcs and the critical temperature at zero magnetic field increase with 
the impurities diffusion potential. The effective interlayer coupling decreases in the presence 
of impurities then the usual superconductivity in the S layer is enhanced. 
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V. EFFECT OF THE IMPURITIES ON THE S-S BILAYER 



In this section, we study the S-S bilayer considering the FFLO modulation and the 
impurities. As predicted irP^ for ferromagnet superconductor multilayered systems, a vr 
state may appear in S-S bilayer under magnetic field. Using the same model as developed 
in the section II, the S-S bilayer is described by the following equations 
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(21) 



where is the superconducting gap in the S layer indexed j 



In the TT— phase, Aq = —A*, the solution of the system (21) is 

1 



/o,o 



1 



and the averaged solution on the angle 



(ni+f72)Ag 



2niQ2(i-^{!^2+f^i)) 

In the clean limit (r — ¥ oo), at zero temperature, the vr superconducting phase appears 
above the critical magnetic field hiow = t — Ag/St and below h^p = t + Ag/St in the limit 
t ^ Aq. As predicted irP^', the modulated FFLO state appears at low temperature and 
maximize the critical magnetic field. Hence, with the FFLO state, the critical magnetic 
fields are hiow,up = ^ =F ^o/^^ in the limit t ^ Aq. The re-entrance superconducting phase 
is enhanced at low temperature by FFLO modulations. The presence of impurities in the 
system may destroy the FFLO state and the re-entrance phase. The FFLO transition should 
meet quickly the first order transition line. Consequently, we will study the infiuence of the 
impurities in the homogeneous case where g = 0. 
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FIG. 12: {h/Tco,t/Tco) diagram for the S-S bilayer in the clean case ((1/2t) /Tcq = 0) (sohd hne) 
and with respectively (l/2r) /Tco = 0.15 (0.194, 0.1, 0.5, 1) in dashed (respectively dotted,dashed- 
dotted and dashed-dotted-dotted) line. The lines h = Aq/2, h = and h = Aq present 

respectively the critical magnetic field for a second order superconducting phase transition, the 
critical magnetic field for a first order phase transition and the FFLO paramagnetic limit for a 
single S layer. The close-dashed line is the FFLO paramgnetic limit in the S-S bilayer in the clean 
limit. 

A. {h, t) phase diagram 



At T = QK in the vr-state without FFLO modulation, the self consistency equation (12) 
becomes : 



The solutions of ( 22 ) are /i*™"^ 




(22) 



up,low 



where h^^j^ow critical magnetic field of the S-S bilayer in the presence of impurities 

(see figure [12] ). The field induced superconducting state is destroyed in presence of 
impurities and cannot be observed if h^^^ = /i™^. We define a critical impurity diffusion 
time Tc = 1/2 2t'^ — a/ 4i(:^ — below which the re-entrance phase totally disappears. 
In the case where t = 2Tco and ho = 0.882Tco then (^)^ ~ 0.194rco. 

The critical magnetic field in the presence of FFLO modulation is plotted in the figure 



19 



12 n the clean limit. The critical magnetic field in the presence of FFLO modulations is the 
upper limit of the critical magnetic field. 

We can see that h^p^^ , the upper critical field in presence of FFLO modulations, cross 
the line h = Aq for t ~ 1.25Tco. It means that the usual superconducting phase is deformed 
only for t > 1.25Tco at T = OK. Then the field induced superconducting phase becomes 
observable. The field induced superconducting phase become totally observable when hf^^^'-^, 
the lower critical field in presence of FFLO modulations, cross the line /i = Aq for t ~ 2.1Tco. 

In the uniform case, we would have to consider the first order transition line. For the S-S 
bilayer, the first order transition line is between the second order and the FFLO transition 
line. The reentrance phase would appear when /i^p, the upper critical field for a first order 
phase transition is above = Ao/\/2 and would be distinguishable if the lower critical field 
in the case of first order phase transition is higher than h. 

B. {h,T) phase diagram 

In the vr-state, the Cooper pairs are formed by two electrons in the different layer. The 
standard superconducting state is only due to the 0— phase and then is not infiuenced by 
the impurities as predicted by the Anderson theorem. The lower and upper critical lines 
merge at field h = t and temperature Tm = vre~'^T|g/ (4t) in the limit t ^ T^q. The field 
induced vr superconductivity is confined to temperature lower than Tm- 

On the phase diagram, we see that the reentrance decreases as the impurity self energy 



is increasing (see figure 13). The reentrance phase totally disappear for 0.194Tco in 

the case where t = 2Tcq. The existence of first order transition line in the field induced 
phase transition could infiuence these results. h{^piow are higher(smaller) than hup^iow Con- 
sequently, the critical impurity diffusion time should be higher than in the case of a second 
order transition. 



VI. CONCLUSION 

To conclude, the proximity effect plays a crucial role in the S-N and S-S bilayers. The 
superconducting critical temperature and the critical magnetic field at zero temperature in 
the S-N and the S-S bilayers depends directly on the interlayer coupling. We demonstrated 
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FIG. 13: {h/Tco, Tc/Tco)phase transition diagram for the S-S bilayer calculated for t = 2Tco with the 
second order transition line in the clean case ((l/2r) /Tco = 0) (solid line) and with respectively 
(l/2r)/rco = 0.01(0.05, 0.1, 0.125, 0.15, 0.175) in dotted (respectively dashed, dashed-dotted, 
dashed-dotted-dotted, close dashed, close dotted) line. The inset presents a zoom of the supercon- 
ducting re-entrance phase around /i 2± f 2± 2T^. 

that at low temperature, a magnetic field induced superconducting phase appears at high 
in-plane magnetic field in S-N bilayers. This field induced phase is originated from the 
compensation of Zeeman effect energy splitting by the energy splitting between the bonding 
and antibonding state electronic levels. This reentrance phase provides the possibility to 
overcome the classical paramagnetic limit and the results of our work give the hints for 
engineering layered superconducting material with very high critical fields. 

In S-S and S-N bilayers, the presence of impurities make the superconducting field induced 
phase more difficult to observe. The impurities produce a broadening of the different energy 
levels over an energy range 1/r which prevents exact compensation. It is possible to define 
a critical mean free pass time over which the re-entrance phase cannot survive. In the S-N 
and S-S bilayer, the critical mean free pass time Tc only depends on the interlayer coupling. 
In S-S bilayer, in the case t ~ Aq then r^T^ ~ 0.25Ao above which there is no possibifity 
to observe field induced phase. Prom an experimental point of view, it could be possible 
with the molecular beam epitaxy techniques to provide a sufficiently large mean free path 
to realize the condition of field-induced phase observation. 

Although we have only treated the Zeeman effect as cooper pair breaking effect, we have to 
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discuss on the orbital pair breaking effect. In the case of multilayered system under in-plane 
magnetic field, the condition for neglecting the orbital effect is given by tH^Qcl/^o < Aq, 
where is the in-plane coherence length and $o = h/2e the superconducting quantum of 
magnetic flux. In the case t ~ Aqwc obtain that H must be lower than Horb — ^0/ i^od). The 
typical values d ~ lOA , ^0 — WOA the corresponding field is extremely large Horb — 200T 
and not restrictive at all as the maximal currently attainable permanent magnetic field are 
GOT. The orbital effect becomes important in layered system in the case t ^ Aq when the 
Pauli limit may be exceeded many times. However ii^, it was demonstrated that the orbital 
pair breaking in layered superconductors are switched off in the high field regime and the 
superconductivity is restored. We may expect that similar situation should be realized in 
S-N and S-S bilayer. 
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